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CN ■ Abstract 

° ■ 

We address the vDVZ discontinuity of the 5D DGP model which consists of a 3-brane residing 



in a flat, infinite-volume bulk. Following a suggestion by Gabadadze [hep-th/0403161] , we 



o 1 implement a constrained perturbative expansion parametrized by brane gauge parameters. We 
^ ■ explore the parameter space and show that the DGP solution exhibiting the vDVZ discontinuity 



corresponds to a set of measure zero. 
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The weakness of the gravitational force has been successfully explained by postulating the 
existence of extra dimensions [1] . The effect of the extra dimensions is a high-energy modification 
of Newton's Law of gravity due to the tower of Kaluza-Klein modes. When the extra dimen- 
sions are of infinite volume, light Kaluza-Klein modes may dominate even at low energies [2-4], 
therefore offering an attractive alternative to dark energy for solving the cosmological constant 
problem. Thus, unlike with finite-volume extra space, 4D gravity is modified at astronomically 
large distances [5]. 4 

The DGP model of a 3-brane residing in a 5D bulk of vanishing cosmological constant [2] is a 
ghost-free, general covariant theory where the 5D graviton mimics a 4D massive graviton on the 
brane. The model appears to be plagued by a van Dam-Veltman-Zakharov (vDVZ) discontinuity 
[7,8] and has attracted much attention [10-12]. In refs.fll], solutions were found interpolating 
between regimes far from and near the Schwarzschild radius by keeping higher-order terms in 
the perturbative expansion. It was thus shown that in the decoupling limit, one recovers the 
standard four-dimensional, weak-field Schwarzschild metric. 

As has been recently argued in [12, 13] for the specific case of D = 5, the breakdown of the 
perturbative expansion at linear order is an artifact of the weak-field expansion itself and can 
be healed by adopting a constrained perturbative expansion. Thus, instead of the incorporation 
of higher-order terms into the linearized treatment, the theory is regulated by a modification of 
the linearized theory itself. After fixing the gauge in the bulk, a residual four-dimensional gauge 
invariance remains on the brane. The graviton propagator is then rendered invertible by the 
addition of a term in the action which would amount to a gauge-fixing term in four- dimensional 
gravity. 

Here, we present a generalized procedure of [12]. We introduce a two-parameter family of 

gauge-fixing terms on the brane. In the decoupling limit they amount to ordinary gauge-fixing 

terms and no physical quantities depend on these gauge parameters. We also fix the gauge in the 

bulk in terms of arbitrary gauge parameters. In the absence of a brane, five- dimensional gauge- 

invariance guarantees that no physical quantities will depend on the bulk gauge parameters. 

We then explore the physical effects of these parameters away from the two extremal limits 
4 See [6] for a slightly different treatment which yields a 4D tensor structure at astronomical distances. 
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(decoupling and absence of a brane). We find that the graviton propagator in general has a well- 
defined decoupling limit implying the absence of a vDVZ discontinuity. The graviton propagator 
exhibits the expected crossover behavior and is found to be free of tachyonic asymptotic states. 
The DGP solution [2] corresponds to a set of measure zero in our parameter space. 

The DGP model describes a 3-brane on the boundary of a five-dimensional bulk-space E. 
The action is 

S DGP = M 3 I d 4 xdy V^G ft (5) + M 2 I d 4 x^ ft (4) (1) 

JS J BY, 

where TZ^ 5 \ (7Z^) is the five- (four-) dimensional Ricci scalar. We adopt the standard conventions 

tjab = diag[H ] ; A, B = 0, . . . , 3, y ; fiv = 0,...,3; i, j = 1, 2, 3. 

Upon varying (1), one arrives at the DGP field equations, which are 

M*Gf B + M 2 G^ B 5(y) = T^ A 5 B 5{y) (2) 

with the linearized solution given by 



M (p 2 + 2m b p) I 3 V 2m bP 

K = -T—e-™^ (3) 
bM rribp 

which is written in terms of the 4D Euclidean momentum and graviton mass 

2 u 2,22,2 

p = -p^p^ = -p Q + Pi =p 4 + Pi 

m b = M 3 /M 2 (4) 

The solution bares a striking resemblance to that of massive gravity where the factor of 1/3 
instead of the Einstein factor of 1/2 signals the existence of a vDVZ discontinuity. In the 
decoupling limit (m& — > 0), 4D Einstein gravity is not recovered and we do not obtain sensible 
dynamics for the longitudinal term with tensor structure of the form p^p v . Although the p^Pu 
term does not contribute at linear level, it does enter nonlinear diagrams. 
Generalizing [12], we define a Constrained DGP Action of the form 

ScDGP = S DGP + S ( V + (5) 

where Sdgp is the DGP action given by (1) and and are gauge-fixing terms in the 
decoupling limit (mj — > 0) and absence of brane (mj — > oo), respectively. Away from these 



two limits (M,M 0), these additional terms no longer simply fix the gauge; they alter the 
boundary conditions. 

We start in the bulk by defining as follows 



= M 3 J d 4 xdy 



2 7 2a 



(6) 



with 



B, = d,hl + ad,h a a -bd a h a , 

By = Ph m (7) 

where a, 7, a, b are arbitrary parameters on which no bulk physical quantities should depend. In 
the absence of the brane, eq. (6) amounts to standard gauge-fixing conditions. In general, the 
a, 7 — ■> limit should be taken at the end of the calculation to ensure that 

-> 

By - (8) 

Next, we define the gauge-fixing term on the brane. For a brane of finite thickness, additional 
terms can arise on the brane world-volume and can survive in the limit of the brane thickness 
tending to zero. In addition, we note that the boundary equations receive no contribution from 
eq. (6) and are invariant under the 4D transformations [12] 

fy^|y=0 — > + d^( u + d u (fi\y=0 (9) 

indicating a residual gauge freedom. With the above in mind, we choose an additional brane 
action contribution 

S-W = \M 2 [ d 4 x^ B V B V (10) 

where 

B v = d^h^ + id u K (11) 

and we assume A > 0. These additional action contributions modify the DGP model by explicitly 
breaking the 4D and 5D coordinate invariance. Adopting this modified DGP model, we next 
obtain and solve the field equations. Varying (5), expanding around a flat background, and 



Fourier transforming, the first-order Einstein equations are as follows. In the bulk, the transverse 
component is 

K - ^f~Kp -\{h y y + «K - b^-hj) = o (12) 

The mixed components are 

id y (p a h a ^ - pji") + p 2 (h m - ^Pu,p a K y ^j + ^p^p a h ay = (13) 
and the components parallel to the brane are 

G$ = (14) 

where 

eg = (p 2 - d y dy)(h^ - %u h a a ) - p^h av - Pl/P a h a , + Vlx p v (K + K) - %Ap 2 H - P a P r K P ) 

+ - -b{p il p u h y y + ap^p u h^ - bp u p a ~h afM ) + ar]^ v (p 2 hl + ap 2 h^ - bp a p /3 h a(3 ) 

+ id v {p v h iiy + p~h vy - 2r] liu p a h ay ) (15) 

From eqs. (12) and (13), we obtain 

~ p f ~ p a p^ ~ \ 

K = b^h a p-ah« + a(h«-^h a ^ (16) 
Plugging these expressions into (15) and assuming the solution is of the form 

h AB (p,y) = h AB (p)e- p y (17) 

we may write Q$ entirely in terms of the AD metric perturbations. Dotting with the momentum, 
we obtain 

pW)g = (1 + a - b)p\p 2 K - p a pP~Kp) (18) 
implying the constraint on the parameters 

l + a-b = (19) 

The vanishing of the divergence, p u G$ = 0, then implies 

p V =p^—^-K(3 (20) 



This is not an additional constraint on the metric. On general grounds, one may argue that 
p v h (W oc p M , hence (20). Using these results, we arrive at the expression 

gj$ = (a - 2 7 - 2a)( PflPu - Vfiu p 2 ) (h a a - ^-hj^ (21) 

leading to a second constraint on the parameters, 

a - 2 7 - 2a = (22) 

At the boundary, the Israel junction condition at y = yields 



M 2 G$ = % v (23) 



where 



G$ = (p 2 + 2m b p)(h iiv - V^K) - (1 - \)(PvP a hau + PuP a h a ») 



+ (1 + 2A0 (p^PuK + y]^p a p r Kp) + 2\e^ v p 2 ~K 

+ 2 imb p - E^j (h° a - ^h a ^j (24) 

Eq. (23) can be solved for arbitrary parameters A,£ and 7. We obtain on the brane 

= — 2 {f, u - L, v Ci + ^C 2 ) f) (25) 

M\p 2 + 2m bP ) I V P 2 ) J 



where 



2mg + 2A(1 + 0(1 - i + 27(1 + i))m b p + A(l + f) Y 

(-1 - 
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6m2 + 4A(1 + 0(1 - 2£ + 3 7 (1 + £))m b p + 2A(1 + V 
(1 - 2A(1 + Q(l + 2 7 (1 + Q)m bP - A(l + Q(l + 2Qp 2 
6m 2 + 4A(1 + 0(1 -2£ + 3 7 (1 + 0)m bP + 2A(1 + V 



(26) 



Notice that the 4.D metric perturbations, when convoluted with a conserved tensor f'^ : 



h^f"" = =2 \% v f'» v - C{Tf '} (27) 

M (p 2 + 2m b p) L J 

are still dependent on the parameters A, £ and 7. Examining the 4D momentum dependence of 

the metric perturbations, we find in the large momentum regime (p ^> m b ) , 

h l ji"< JV ~ {T^T'^ - iTT'j (28) 

M p 2 1 * 



recovering 4D Einstein gravity, and in the small momentum limit (p <C m b ), 

* ^ " 1^'} (29) 

exhibiting 5D behavior, as expected. Notice that in both limits, the transverse components of 
the metric on the brane are independent of the parameters A, £ and 7. In the intermediate range, 
the propagator smoothly switches from the 4D expression (28) to the 5D expression (29) as the 
momentum decreases. This crossover behavior depends on the parameters A,£ and 7. 

In the decoupling limit, m b — > 0, the graviton propagator yields the standard 4D Einstein 
solution on the brane demonstrating the absence of a vDVZ discontinuity. This is the case in 
the entire parameter space except for a set of measure zero defined by 

£ = -1 (30) 

For this special choice, the parameters become true gauge parameters throughout the entire 
range of momenta. We obtain 



M (p 2 + 2m b p) I 3 V 2m 6 p 

K(p,y) = — re-" 

6M m b p 

h m (p,y) = (31) 

which is independent of a, 7. Also, the constraints = B 5 = for general a, 7 showing that 
they represent gauge-fixing conditions. This is the standard DGP model [2]. 

It should also be noted that for the choice of parameters A = 1, £ = —1/2, we recover the 
model proposed by Gabadadze [12], 



M (p 2 + 2m b p) I 2 (p 2 + Qm b p)~ 

Kv = -^(~K-^hap) (32) 



P V P 

in the a, 7 — * limit. 
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We next wish to examine the poles of the propagator. Taking the 7^0 limit, the transverse 
part of the propagator (27) can be written in a form explicitly revealing its pole structure, 

VT V = =2 T^f > - C{p)ff (33) 

M (p 2 + 2m b p) 3M 



where 



C(P) = ., , + 



p 2 + 2m^ 2(c + — c_ 



(p 2 + c + m b p) (p 2 + c-m b p) 



The location of the poles is determined by the coefficients 

1 



c± = c±(£, A) = 



1 + e 



1 -2f.± \/(l-20 2 -| 



(34) 



(35) 



For p 3> m&, C(p) w ^ an d we recover the 4D expression (28). The poles are significant 
for momenta p < m^. As was shown in [12], the p = —2m b pole lies on the second Riemann 
sheet in the Minkowski four-momentum complex plane, where p 2 = s exp(— in), s = p^f. 
This pole corresponds to a non-physical resonance and indicates an intermediate, metastable 
state. This can be seen from the p = ±^/— s dependence of the propagator which indicates that 
the propagator is multi-valued and the complex s-plane has two sheets with a branch cut on 
the positive real axis. For the choice of p = \/—s, we obtain a non-physical resonance and a 
propagator which decays with the bulk coordinate. 

The other two poles are located at p — —c±m b and depend on the parameters £ and A. In 
the (£, A)-plane, above the curve 

A = ^— - (36) 

(l-2£) 2 1 ; 

both poles lie on the negative real axis in the complex s-plane, since c± G M. Moreover, c± > 
for — 1 < £ < 1/2. In this strip, the two poles are in the second Riemann sheet (corresponding 
to the choice p = t/—s) and are thus unphysical. In the special case £ = —1/2, A = 1, the pole 
at p — —C-rrib coincides with the pole at p — —2mi,; this is the Gabadadze model [12]. As we 
approach the curve (36), the two poles come together. Below the curve (36), c± become complex 
and c + = c*_. In this case the poles are no longer on the real axis; we obtain a resonance with a 
momentum independent decay width, in addition to the pole at p — —2m b . 

As £ — > —1, the two poles p = c±m b become infinite and C{p) — > (p 2 + 2m b p)~ x . This is 
a singular case; dependence on the A parameter disappears and the propagator turns into the 
DGP expression (31) [2] which is plagued by the vDVZ discontinuity. 

8 



To the left of £ = — 1 (as well as for £ > 1/2), both c± < 0; therefore, the poles p = —c±rrib 
are tachyons, signaling instability of the solution. Were we to choose p = — s, instead, we 
would place these two poles on the second Riemann sheet, but then the third pole at p — —2rrib 
would turn into a tachyon. 

The above results are illustrated by the two-dimensional plot of the (£, A) parameter space 
in Figure 1. 

In summary, we generalized the constrained perturbative model of [12] and calculated the 
graviton propagator. The first-order contribution to the perturbative expansion depended ex- 
plicitly on parameters which are gauge parameters in the bulk (in the absence of a brane) and 
on the brane (in the decoupling limit), respectively. These parameters determine the details of 
the distance-dependent, crossover behavior of the propagator and the position of the poles of the 
graviton propagator. At low momenta, we obtained a 5D behavior whereas at high momenta we 
recovered 4D gravity demonstrating the absence of a vDVZ discontinuity. In addition, we found a 
range of parameter values which yielded non-physical resonances corresponding to intermediate, 
metastable states. For a special choice of parameters (representing a set of measure zero in the 
parameter space), we recovered the standard DGP model [2]. This choice represented a set of 
measure zero in the parameter space which is plagued by the vDVZ discontinuity. 

It would be desirable to understand the origin of these parameters better as physical quantities 
depend on them. They may represent different physical setups (embeddings of the brane in the 
bulk) or "schemes" (similar to QCD) which are artifacts of the perturbative expansion and would 
be resolved once higher-order terms are included. We hope to report on progress in this direction 
shortly. 
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Figure 1: The two-dimensional (£, A) parameter space. Above the curve (36), all poles of the 
propagator are real. Within the strip — 1 < £ < 1/2, only unphysical resonances appear; outside, 
we have tachyons (instability). Below the curve (36), we have one real pole and a resonance with 
momentum independent decay width. The DGP model [2] is represented by the line £ = — 1; the 
Gabadadze model [12] by the point £ = —1/2, A = 1. 
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